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FOREWORD 


This final report presents a summary of the analytic 
model developments and results of analyses performed 

in the Structural Dynamic Analysis of the Space Shuttle 
Main Engine. This work was performed by the Structures 
& Mechanics Section of the Lockheed Missiles & Space 
Company, Htmtsville Research & Engineering Center, Hunts- 
ville, Alabama, for the National Aeronautics and Space 
Administration, George C. Marshall Space Flight Center, 
Marshall Space Flight Center, Alabama. The Contracting 
Officer's Technical Representative for this study is 
Mr. Larry A. Kiefling, ED22. 

Thermal and structural models of critical SSME powerhead 
components developed and applied in this study are docu- 
mented and analyses described in this final report. These 
documented models, installed on MSFC computers, provide 
a data base for future SSME engine uprating studies and 
modifications. 
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1. INTRODUCTION 


The objective of this structural dynamic analysis is to support dev- 
elopment of the SSME by evaluating components subjected to critical dynamic 
loads, identify significant parameters and evaluate solution methods. Engine 
operating parameters at both rated and full power levels (RPL and FPL) were 
considered. Detailed structural dynamic analyses of operationally critical 
and life-limited components have supported the assessment of engine design 
modifications and environmental changes. Engine system test results which 
becamct available during the work period have been utilized to verify analytic 
model simulations. 

Tlie SSME Main Chamber Injector Assembly shown in Fig. 1 is an assembly of 
600 injector elements which are called LOX posts. It is this critical compo- 
nent which is the focus of the analysis effort described herein. 

The overall LOX post analysis procedure is shown in Fig. 2 where the 
analysis flow proceeds from thernuil, structural static and dynamic loadings 
to the nonlinear APSA stress analysis. Later, in order to develop a more 
detailed view of the nonlinear cyclic loads, a BOPACE 3-D model was developed. 
The analysis is concluded with the fatigue life assessments. These procedures 
were, of course, often repeated for the evaluation of design and material 
changes. 

A sectional view of the assembled LOX post is shown in Fig. 3 where each 
component is identified and a set of static applied lateral deflections are 
summarized. 
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Fie. 1 - SSME Main Chamber Injector Assembly 
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Fig. 2 - SSME LOX Post Analysis Procedure 
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2. ANALYTIC MODEL DEVELOPMENTS 


The finite element math models de\eloped for the LOX post analysis are 
summarized as follows: 

1. SPAR LOX Post Assembly Model 

• Thermal Version 

• Stress Version 

• Dynamic Version 

2. SPAR Post/Retainer Detail Model 

• Thermal Version 

• Stress Version 

3. APSA Post/Retainer Detail Model 

• Stress Version (Nonlinear) 

4. BOPACE 3-D Post/Retainer Detail Model 

• 60 Deg Symmetric Segment 

5. BOPACE 2-D Inertia Weld Region Model 

• Stress Version 

6. BOPACE 3-D Inertia Weld Region Model 

• Thermal Version 

• Stress Version 

7. SPAR Inertia Weld Region Model 

• Thermal Version 

• Stress Version 
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Two different post configurations were analyzed: (1) a one-piece 316L 

stainless steel post, and (2) an improved two-piece post with a Haynes 188 
tip. Modified versions of the LOX post models were used to simulate posts 
from Rows 13, 12, 11 and 10; however, most analyses were concentrated on the 
Row 13 "Baseline" configuration which was investigated both with and without 
the outer flow shield. 

2.1 SPAR LOX POST ASSEMBLY MODEL 

The SPAR LOX Post Assembly Model Components are shown in Fig. 4. A lay- 
out drawing of the complete SPAR LOX Post Assembly Model and its 46 cross sec- 
tions with numbered node locations are presented in Appendix A. Also shown 
is the one-half symmetric 13th row flow shield model details. 

Assembled and exploded views of the LOX Post Assembly Model Components 
are presented in Appendix B. 

This SPAR assembly model of the 13th row LOX Post contains 3388 nodes 
(2964 in the post assembly and 424 within the flow shield). The structural 
model possesses 10164 dej^rees of freedom and contains 2400 three-dimensional solid 
elements. The same element mesh is used for the thermal and structural models. 

The Row 12 model is developed from the baseline by deleting the shield and 
adjusting the distance between the upper post inertia weld joint and the sec- 
ondary plate. 
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2.2 SPAR POST/RETAINER DETAILED MODEL 

The SPAR Post/Retainer Detail Model was developed to investigate both 
thermal and stress levels in the post-to-retainer thread area where fatigue 
cracks had been observed after test firings at both RPL and FPL levels. The 
Post/Retainer Thermal/Structural Model is shown in Fig. 5. Additional views 
of the model and its components are presented in Appendix C. The detailed 
model consists of 2040 nodes. 

2.3 APGA POST/RETAINER DETAIL MODEL 

The APSA Post /Retainer Detailed Model was developed to Investigate the 
nonlinear effects encountered in the post-to-retainer threads which resulted 
in cracks forming in the threads and fillet region of the post. The APSA 
model shares exactly identical boundary nodes with the SPAR assembly model so 
that computed temperature prof.les from the SPAR thermal analysis could be 
mapped directly onto the APSA detailed model for stress analyses. The APSA 
axisymmetric model is shown in Fig. 6. Because of the detail involved, only 
the first three threads are modeled. The APSA model has 1500 nodes, 2958 degrees 
of freedom and is constructed using 1359 axisymmetric elements. An evaluation 
of the stress variations due to load sharing between the threads indicated 
that three threads represented an adequate modeling approach. An enlarged 
view of the threaded area is shown in Fig. 7 where thread contact stresses 
are simulated through the use of "slip" elements at the interface. 


2-4 


LOCKHEtO- HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D734A90 



2-5 


IcVKutUt MUMSVilU HEStARCM & INGlNtfHlNG CENUR 


Hg* 5 - SPAR Post/Retainer Detailed Model 


LMSC-HREC TR D784490 


C 0000X10 


O S. 0000X10 


•• 0000X10 


. 3 0000X10 


? 0000X10 


I 0000X10 



/ OlHfNSIOS INCHES 


Fig. 6 - APSA LOX POST/Retainer Detailed Thread Model 
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Fig. 7 - APSA Post/Retainer Detailed Model - Thread Region 
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2.4 BORACK )-0 I'OST/RRTAINKR nKTAIl.KD MOPEL 

The BOl’ACK J-P Post /Retainor Potaile-i Moilol is maJo up of a hO-det; synimo- 
trio sonnaMU of iho lliroail rojtion with 810 nodes and 405 eight-node brick ele- 
ments. Elastic-plastic properties are included for JlbL stainless steel and 
Haynes 188. At the 0-deg axial plane nodes are numbered 1 through 270, at the 
JO-deg plane, nodes are numbered 1001 through 1270 and at the bO-deg plane, 

2001 through 2270. The current run stream is configured for processing on the 
IBM ibO computer system at MSEC. I’he model comprises the same post/retainer 
thread reglmi used in the SPAR and APSA detailed models. This BOPACE simula- 
tion, however, permits multiple elastic-plastic load cycling not available in 
the other models. Petailed model plots have been generated and are available. 

2.5 BOPACE 2-P INERT 1 \ WEl.P REOION MOPKL 

To evaluate the stress concentrat ion effects which exist at the l.OX post 
inertia weld Joint, a refined local two-dimensional model was constructed. This 
2-P Inertia Weld Region Model is presented in Fig. 8. The refined region at 
the outer weld plane between the 31bl. and Haynes 188 materials permits accurate 
evaluation of the stress concent rat ion due to the chemical etching performed 
after welding. An example of the method of simulating a variable depth of 
etching is slu^wn in Fig. ‘J where the appropriate elements are removed to simu- 
late a O.OOb-in. etch depth. I'he refined area around the weld is nuide up of 
ten lavers each 0.001 it\. thick. Both uniaxial and bending loadings were 
applied for the evaluation ot the stress concentration as a function of etch 
depth . 
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Fig. 8 - BOPACE 2-D Weld Model (Unetched) 
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Fig. 9 - Model Loading and Element Layout for 0,006 in. Etch Depth 
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2.6 BOP ACE 3-D INERTIA WELD REGION MODEL 

Figure 10 presents the detailed BOPACE 3-D model of the inertia weld 
region built to investigate the post-yield performance of the weld joint. 
This 0.115- in. long, 5-deg segment model is built of eight-node brick ele- 
ments. Tect data reveal significant deformation and cracking near the weld. 
Thermal and structural load conditions were generated for this model. 
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Fig. 10 - BOPACE Model of the LOX Post Inertia Weld 
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2.7 SPAR INERTIA WELD REGION MODEL 

In order to perform the required thermal analysis to define accurate weld 
region environments a detailed SPAR model of the upper post and inertia weld 
region was developed. The upper LOX post assembly and its components are 
shown in Fig. 11. The detailed SPAR model of this region is shown in Fig. 12. 
This 360”deg three-dimensional model contains 1932 nodes. All SPAR models 
are configured for processing on the MSFC Univac 1100 computer system. 
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3. LOX POST ANALYSIS DISCUSSION 

This discussion is subdivided into three areas which are arranged as 
follows: 

1. LOX Post Assembly Analysis 

• Thermal 

• Structural 

• Dynamic 

• Fatigue 

2. LOX Post Inertia Weld Region Analysis 

• Thermal (at FPL) 

3. Inertia Weld Region Stress Concentration 

• Uniaxial Loading 

• Bending Loading 

3.1 LOX POST ASSEMBLY ANALYSIS 

The RPL thermal analysis was performed on the SPAR main chamber LOX post 
assembly model using the environmental temperatures and film coefficients 
identified in Fig. 13. The analysis results from the assembly model thermal 
run were applied to the Post/Retainer Detailed Model which was used to dev- 
elop the temperature distributions summarized in Fig. 14 for the one-piece 
316L stainless steel configuration and in Fig. 15 for the Haynes 188 two- 
piece LOX post. The radial temperature results are presented at five cross- 
sections through the post/retainer thread region. 


I 
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Fig. 13 - LOX Post Thermal Load Input*- for All Models 












Fig. 14 - Radial Temperature Variation in Post/Retainer Thread Area for 316L 
One-Piece LOX Post 
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Fig. 15 - Radial Temperature Variation in Post-Retainer Thread Area for Haynes 188 
Tip LOX Post 
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The exact RPL temperature profiles from the SPAR Post/Retainer Detailed 
Model were nuipped onto the APFA Thread Model as shown in Figs. 16, 17 and ]8. 

In addition, the SPAR model was used to develop the displacements due to non- 
axisymmetric loads (i.e., pressure and transverse motion) that exist at the 
detailed model boundary nodes. These loads are summarized for two cases in 
Fig. 19 for the one-piece 316 L configuration and in Fig. "lO for the two- 
piece Haynes 188 LOX Post. 

Vibration induced loads were calculated using the SPAR LOX Post Assembly 
Model. The general procedures used are outlined in Fig. 21, where the SPAR 
model response is forced to reach measured levels of strain in order to access 
the magnitude of the vibration- induced loads. 

Figure 22 summarizes the first three norrial modes of each configuration 
both with and without the flow shield. The resultant dynamic displacements 
were determined and applied to the APSA model as shown in Fig. 23 for the one- 
piece 316 L LOX Post and in Fig. 24 for the two-piece Hynes 188 version. 

All these corabxned loadings were evaluated in the stress analysis of the 
APSA mode]. Peak resultant stresses were identified at points D and C in Fig. 

25. These maximum stresses and strains are summarized in the table in Fig. 

26, where the peak effective strain of 3.2% was experienced in the 316 L 
stainless steel post thread area. 

Stress-strain diagrams of the first one-half cycle of combined static, 
thermal and dynamic loading effects for each of the three peak value cases 
are presented in Figs. 27, 28 and 29. Here, it can be seen that for all cases 
the state of stress is well within the plastic or post-yield range for each 
case. Inspection of the S/N curves for the two materials reveals the following 
endurance limits: 

316 L (Fig. 30) 0 = 16 ksi 
Haynes 188 (Fig. 31) c = 36 ksi 
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Using Goodman's correction for an equivalent stress, neither of the two 
configurations reaches its endurance limit; 

316 L: 41% of endurance limit (See Fig. 32) 

Haynes 188: 20% of endurance limit (See Fig. 32) 

It should be noted however that the one-piece (316 L) peak stresses are 
experienced in a tension field while tnose of the two-piece (H-188) post are 
in a compression stress region. The Haynes 188 nuiterial in the two-piece 
design can be concluded to be a superior configuration from this fatigue 
analysis. 
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Fig. 16 - APSA Thread Fillet Model Temperature Distribution 
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Fig. 17 - APSA Model Temperature Contours 
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Fig. 18 - APSA Model Temperature Contours 
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Fig. 19 - One-Piece 316 L LOX Post Displacements Due to Non-Axi symmetric Loads 

















Fig. 20 - Two-Piece (Haynes Tip) LOX Post Displacements due to Non -Ax i symmetric Loads 
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Piece 316 L LOX Post Displacements due to Measured Vibration Loading 
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Fig. 24 - Two-Piece (Haynes TT^.; LOX Post Displacemnts Due to Measured Vibration Loading 





















APSA THREAD MODEL - DETAILED STRESS POINTS 




MAXIMUM STRESS AND STRAIN IN ROOT AREA OF THREAD ON LOX POST 

(Combined Thermal, Static and Dynamic Displacement Loads) 
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Fig, 26 - APSA Thread Model Results 















Post Material = 316 L 
Evaluation Element 461 
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Fig, 27 - Stress-Strain Diagram for Combined Static, Thermal and Dynamic Loads 
at RPL Conditions 






Post Material = Haynes 
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Fig. 28 - Stress-Strain Diagram for Combined Static, Thermal and Dynamic Loads 
at RPL Conditions 








Pobt Material = Haynes 188 
Evaluation Element 499 
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Fig. 29 - Stress-Strain Diagram for Combined Static, Thermal and Dynamic Loads 
at RPL Conditions 
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Fig. 31 - Fatigue Life Curve for Haynes 188 
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Fig. 32 - Calculation of Effective Dynamic Stresses 
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3.2 LOX POST INERTIA WELD REGION ANALYSIS 

Using the existing SPAR detailed upper post model, a complete FPL thermal 
analysis was performed using environments and film coefficients specified by 
MSEC. A view of the SPAR model in the weld region is shown in Fig. 33 for 
both the undeformed and deformed conditions. The complete distribution of 
resulting FPL temperatures at the weld region are shown in Fig. 34. Signifi- 
cant temperature gradients can be seen between the post ID and OD surfaces. 

Tlie two different materials in the area are identified. Appropriate tempera- 
ture dependent conduction properties were used for each in the analysis. 

Using the thermal environments determined in the above analysis, a detailed 
BOPACE elastic-plastic stress analysis was performed with the 5-deg BOPACE 
Weld Region Model. The resultant strain levels iii each material are shown 
in Table 1 at the weld line. The effective strain ranges due to heating 
from room temperature to FPL and then back again are shown in the table. These 
values v'ere used to assess the low cycle fatigue life of both the 316 L and 
Haynes 188 LOX posts. 

The 316 L low cycle fatigue life at 0.32% strain range is seen in Fig. 35 
to be approxinuitely 10,000 cycles. Similarly, the Haynes 188 at 0.53% strain 
range has a low cycle life of approximately 2000 cycles as shown in Fig. 36. 

Thus, in summary, the predicted FPL cycle does produce thermally induced 
plastic strain in the inertia weld region of the LOX post for both materials. 

The thermal strain cycle is more severe in the Haynes 188 post than in 316 L 
due to the lower icsistance to low cycle fatigue of the Haynes material (2000 
cycles in Haynes versis 10,000 for 310 l). However, in either case, adequate 
low-cycle fatigue life levels are predicted. 
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Table 1 

BOP ACE INERTIA WELD MODEL 
(COMPUTED THERMAL STRAINS - FPL CYCLE) 


LOX Pest 
Material 

1 

Effective Strain, Outer Surface (%) 

Effective Strain 
Range (%) 

at FPL 

Residual at R.T. 

316 L 

+0.48 

+0.16 

+0.32 

(As Welded) 




Haynes 188 

+0.69 

+0.16 

+0.53 

(As Welded) 
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Pig. 35 - Low Cycle Fatigue Curves for 3j.b L 
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Fig. J6 - Low Cycle Fatigue Curves for Haynes 188 
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3.3 INERTIA WELD REGION STRESS CONCENTRATION 

An analysis investigation was performed of the stress concentration 
factors which exist in the upper post inertia weld region of the main stage 
LOX post after etching is performed. A two-dimensional BOPACE analytic 
model of the weld region was developed as described in Section 2.3. Here the 
two material regions are divided by the diagonal weld line representing one 
side of the tubular post. The most likely crack forming location is the out- 
side region of the weld which has been refined to simulate ten layers (each 
0.001 in. thick). 

A series of load cases was analyzed using uniform stress loads. Tension 
loads were applied to the left boundary of the model while the right side was 
held stationary. Etch depth variations from zero to 0.009 in. were evaluated 
by the sequential removal of element levels as shown in Fig. 37 where the 
0.006 in. depth is depicted. 

Stress resultants were determined at Points 1 and 2 as shown in the enlarged 
model view in Fig. 37. Stresses for each etch depth were compared with the 
uniform stress field of the unetched mouel. Stress concentration factors, K^, 
were determined at each location for each etch depth. Results for the elastic 
range of factors are presented in Fig. 38 as a function of etch depth where a 
factor of 1.0 represents the zero or unetched model. 

Similarly resultant elastic-plastic stress concentration factors are pre- 
sented versus percent load for various etch depth in Fig. 39 where 100% corres- 
ponds to a uniform stress equal to the 316 L material yield strength for the 
unetched model. It can be seen that as etch depth Increases the weld region 
stress concentration causes plastic material performance at earlier load levels. 
Unfortunately, this lowering of the local stresses is accompanied by work 
hardening of the material which in turn raises the yield stress and causes the 
general maintaining of elastic ILp values. For this reason, it is recommended 
that application of values be made at the higher, l.e. elastic, values and 
not the lower plastic levels. Advantages of lower etch depths are clearly 
demonstrated in these results. 
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An additional evaluation was performed loading the BOPACE model as 
shown in Fig. 40. This load distribution was derived from a pure bending 
moment load on the LOX post tubular crossection. Ihe 100% load level was 
adjusted to produce 316 L yield stresses at the etch gap location in the 
unetched model. Resultant values versus depth of etch at Points 1 and 
2 are shown in Fig. 41. The elastic-plastic performance is presented in 
Fig. 42 versus percent load. Peak values for this bending moment case 
approach a value of 2.5; whereas for the earlier pure tension case values 
approach a peak of 2.0. 
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Fig. 37 - Inertia Weld Region Enlargement 
Showing Stress Evaluation Points 
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Data Reduced from Elastic Strain Region 
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Fig. 38 - LOX Weld Stress Concentration Factor Versus 
Depth of Etch 
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Fig. 41 - LOX Post Weld Concentration Factor vs Depth of Etch 
for Bending Moment 
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K^, at Etch Depth 



Fig. 42 - LOX Post Weld Area Stress Concentration Factor vs 
Percent Load for Bending Moment 
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Post Inner Surface Convective Elements 
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